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Abstract: Theories beyond the standard model such as string theory motivate low en- 
ergy effective field theories with several scalar fields which are not only coupled through a 
potential but also through their kinetic terms. For such theories we derive the general for- 
mulae for the running of the spectral indices for the adiabatic, isocurvature and correlation 
spectra in the case of two field inflation. We also compute the expected non-Gaussianity 
in such models for specific forms of the potentials. We find that the coupling has little 
impact on the level of non-Gaussianity during inflation. 
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1. Introduction 

String theory and theories with supersymmetry usually contain many scalar fields which 
can play an important role in the early Universe. For example, in string theory they often 
describe the dynamics of extra spatial dimensions and other degrees of freedom living in 
higher dimensions. These scalar degrees of freedom couple to matter fields propagating in 
the three large dimensions we perceive, see for example [l]-[3]. If these extra dimensions 
exist, they will have had an influence on the evolution of the universe at some point. It 
is usually thought that in particular the dynamics of the very early universe is affected by 
the existence of extra dimensions. As such, they will alter the predictions of inflationary 
cosmology. If apart from the inflaton field(s) other scalar fields are present, they will 
generally alter the evolution of the field(s) driving inflation and affect the production of 
cosmological perturbations. 
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Generally speaking, the existence of multiple fields during the inflationary epoch mod- 
ifys some of the single field predictions. For example, apart from the usual adiabatic 
perturbations produced during inflation, there could be isocurvature (entropy) perturba- 
tions produced, whose existence is constrained by observations (see, for example [4]-[13] 
and references therein). 

Related to this, in single field inflation, the curvature perturbation on constant energy 
density hyper surf aces, £, is constant on super-horizon scales and can be evaluated at the 
horizon crossing. However, in the presence of multiple fields during the inflationary epoch, 
( does not remain constant and varies since the non-vanishing isocurvature perturbations 
act as a source term for the change of ( [5, 6]. 

In order to distinguish between inflationary models, cosmologists need to extract as 
much information as possible from the data. Usually the spectral index and the scalar- 
tensor ratio are used to distinguish between inflationary models. With the advent of 
precision data (in particular the data obtained by WMAP[14]), the running of the spectral 
index can be added as an additional quantity in order to distinguish between the models. 
One can expect considerable improvement with future data coming from the PLANCK 
satellite and the mapping of the large scale structures in the universe, as well as a better 
understanding of small scale clustering as obtained from the Lya forest. 

Another potential observable which could distinguish between inflationary models is 
the amount of non-Gaussianity generated during inflation (see [15] and references therein). 
Although it is small for many inflationary models (e.g. [16]- [18]), there are examples in 
which perturbations show a considerable amount of non-Gaussianity [19]-[26]. 

In this paper we study the slow-roll regime of generalised two-field inflation, where the 
scalar fields are also coupled through kinetic terms [5, 6, 9]. In general, a multi- field action 
can be given by 



/ 



d x \J — g 



M^-I^d^W -w{4>) 



(1.1) 



For simplicity we consider a system of two scalar fields, whose dynamics are governed by 
the following action: 



= Jd A x^-g Ml^- l -g^d^d^-\e 2h{lp) g^d, X dvX-W{v,x) 



(1.2) 



In this expression Mp = l/\/8irG is the reduced Planck mass. Note that, in this case, Qu 
is symmetric. 

We obtain the general formulae for the running of the spectral indices (adiabatic, non- 
adiabatic and correlated), following [27] and [28]. The running is affected by the presence of 
the coupling and therefore future experiments constraining the running of the spectral index 
will give vital information about the existence of fields which couple non-trivially to the 
field(s) driving inflation. Furthermore, we derive the general formulae for non-Gaussianity 
and apply them to some specific models. 

The paper is organised as follows: In the next Section we summarise the equations 
governing the background and perturbation evolution. In Section 3 we derive the general 
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expressions for the running of the spectral indices (for the adiabatic and isocurvature 
power spectra as well as the correlation spectrum). In Section 4 we derive the expected 
non-Gaussianity for two cases. In Sections 5 and 6 we then apply our results to specific 
examples. Numerical results are presented in Section 7. Our conclusions can be found in 
Section 8. Useful formulae are collected in the Appendices. 



2. Background and Perturbation Equations 



In this section we set our notation and briefly review some results found in previous work. 
It will closely follow [27] and [28]. The equations of motion for the two scalar fields in a 
Friedmann-Robertson- Walker spacetime follow from the action (1.2) and are given by 



(p + + W v = V 2b x 2 , 
X + (3H + 2b ip <p) X + e~ 2b W x = 0, 



(2.1) 



where b v = 9b g^ , W v = 9W g^,' x ^ , etc. Einstein's equations lead to 



H 



3MJ, 



2Mp 



^ 2 + e 2b x 2 



W 
a 2 



2M 2 p ' 



We define e-folding number, N, as 



N(t e ,Q = ! Hdt. 

J * 



(2.2) 



(2.3) 



It is useful to separate the perturbations into components of adiabatic and isocurvature 
modes. Following [27, 29], we define the average (adiabatic) and orthogonal (entropy) fields 
a and s as 



da = cos 9d(p + sin 9e b dx, 
ds =e b cos 9d\ — sin 9d(p, 



(2.4) 



with 



cos ( 



sin I 



e b X 



vV + e 2b x 2 



These fields satisfy the equations of motion 



(2.5) 



a + 3H& + W a = 0, 



b^a sin 9, 



(2.6) 



where 



W„ 



W 9 cos 9 + e~ b W x sin0 
-W^sm9 + e' b W y cos9 



(2.7) 
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We now turn to cosmological perturbations. We work in the longitudinal gauge, in 
which the perturbed metric has the form 



ds 2 = -(1 + 2<S>)dt 2 + a 2 (l - 2$)dx 2 , 



(2.8) 



where a = a(t) is the scale factor and = 3>(t, x) is the metric perturbation. We now 
calculate the perturbations in the adiabatic and entropy fields. Instead of working with the 
field perturbation 5a, it is more convenient to work with the Sasaki-Mukhanov variable 
Q a = 5a + jj<&. The equation of motion for Q a is given by 



Q a + 3HQ a + 



k 2 



^ + w aa + e—- 



M 2 a 3 



a 3 a 2 
H 



\ a H a 



, W y e- b 

bipip — sin( 

a 



Qa 



(2.9) 



whereas for 5s it is given by 
~k 2 



5s + 3H5s + 
with 



+ W ss + W 2 - b w a 2 + b% git) + b v f(t) 



5s 



k 2 $ W s 
^2^7^' 



(2.10) 



g (t) = -CT 2 (l + 3sin 2 6»), 

f(t) = W v (l + sin 2 9) - 4W S sm6. 



(2-11) 



We have taken the notation of [28] and use 



W aa = W w cos 2 9 + V vy sin 29e~ b + V yy sin 2 9e 



-2b 



'H>X 



XX ' 



W ss = W w sin 2 9 - V m sin 29e~ b + V xx cos 2 9e~ 2b . 



(2.12) 
(2.13) 



We define the curvature and isocurvature fluctuations as 



( = ^-Qa, S = ^-5s, 
a a 



(2.14) 



then the time derivative of ( is related to S by 



Ha 2 a 



(2.15) 



During inflation, the fields are assumed to follow the slow-roll limit, 



— a cos — -g^, 



X = <Jsm6e b ~ -^e 2b , 



h 2 (^x)^^tW(^ x ). 

The slow-roll parameters are defined in Appendix A for convenience. 



(2.16) 
(2.17) 
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In the slow-roll limit and on large scales, we find the evolutions of curvature and 
isocurvature perturbation can be written in terms of slow-roll parameters, 



where 



C ~ BHS, S ~ -jHS, (2.18) 
-2r] as - sign(6 v ,)sign (j^J \/^sin 2 9, 

-i] aa + 2e + rj ss + i sign %) sign (j^f"j \f^x sin 6 cos 9 ( 2 - 19 ) 
+ ^sign(6^)sign (^jfj V^i 1 + sin2 

Generally speaking on large scales the evolutions of adiabatic and isocurvature fluctuations 
follow the following set of equations 



B 
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C = a(t)H(t)S, 
S = 5(t)H(t)S. 



(2.20) 



We use the formalism of transfer matrix [4], 



C(*) = f 1T cs| C(*. 

S(t) ) [OTss { S(U 



(2.21) 



where 



(2.22) 



T SS (U,t) = exp ^ 5{t')H{t')dt' 
T CS (U,t) = f ' a{t')H(t')T S s(U,t')dt' 
Then the power spectra are 

P c = (1 + Tls) Vd = V C U(1 + cot 2 A), (2.23) 

Vs = T 2 ss V c U, (2.24) 
Vc = T cs TssV c U, (2.25) 
where the cross-correlation angle A is 

Vc 

cosA = — ===. (2.26) 

We note that T^g = cot A and A has range < A < it giving positive and negative 
correlation depending on the sign of cos A. The spectral indices are defined as 
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The spectral indices are best written in terms of slow-roll parameters. They have been 
calculated in [28] and are given by 



n ( - 1 = -6e + 4e(cos A) 2 + 2rj a(T (sm A) 2 + 4q as sin A cos A + 2v ss {cos A) 2 



+ sign (b^ ) sign I — j ] ^/e h e v (l + sin 2 9) cos 2 A 



+ 2 sign(b (/3 )sign I —A j ^/ebe x (sin 9) 2 sin A cos A 



— sign(b ((3 )sign ( —A j ^/ebe x sin6> cos sin 2 A, 




(2.28) 



n s - 1 = -2e + 2r/ ss + sign (b^) sign 




^/^(l + sin 2 #), (2.29) 



and 



rtc — I = ns — 1 — B tan A. 



(2.30) 



In these expressions, all the slow-roll parameters are evaluated at horizon crossing. 



3. Running of the Spectral Indices 

To calculate the running of the spectral indices we need to know the second-order slow-roll 
parameters and the derivative of first-order slow-roll parameters and the transfer functions. 
For convenience, we list them extensively in Appendix A. 

When ^ we have 14 parameters at second-order in general, as listed in Table 1. 
However, considering the component fields (cr,s), H, T$s and £ 2 S6 , do not appear in the 
running of spectral indices, which reduces the parameters to 11. Therefore the running 
spectral indices are given by: 6, A and 9 slow-roll parameters. We define the running 
spectral indices ax as 



(III v 
dink 



x = c,s,c. 



(3.1) 
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To leading order in slow roll, the runnings of the spectral indices after inflation are: 



drit 



a^ 



dink \ d In k 



dn, 



C 



+ 



6=0 



16 



5 - 3 cos 26* + 8(1 - cos26>)(sin26>sin2A + cos 26* cos 2A) 



-2 sin 26»(2 - 3 cos 26 + cos 2 29) sin 4A - {3 + cos 2(9(3-6 cos 29 + 2 cos 2 26*)} cos 4A 



1 

8 L 



1 + cos 29 + 2 cos 2 A + ( 1 - cos 29) (sin 29 sin 2 A + cos 29 cos 2 A) 
+ 2 cos 3 6» + cos 6»(3- 2 008 26*) cos 2A + cos 9(— 2 + cos 29) cos 4A 
+ 2sin 3 6>(2 sin 2A — sin 4A) sign(6 (f3 )sign 

+ 



6>e 



K 

V 



e&ee 



a s = 



dn s 
dink 



| - cos0 + sin6»(l - cos 26>)(- sin 2A + ^sin4A) + - cos 0(-3 + 2 cos 26>) cos 2A 
+ ^cos6 , (2 - cos26 ) )cos4a|t? CT(T + | ^ sin 6> + i cos 6>(3 + cos 26>) sin 2A 
- cos 6»(2 - cos 29) sin 4 A - - sin 9(6 + cos 29) cos 2 A + sin 9(1 - cos 29) cos 4a|?7 cts 
+ | - cos 9 + - sin 9(3 + cos 26>) sin 2 A - sin 3 9 sin 4A + cos 3 9 cos 2 A 
+ i cos 6>(-2 + cos 26>) cos 4a|?? ss sign(6 v )sign (yf^J 

I "77^- I + sin 2 6* cos 4 6>e b e - - (1 + sin 2 9) cos 2 6>£ 6 e + [2 cos 9e - ( 1 + sin 2 9) cos % 
\dlnkj b=0 2 

+ sin6 ) (-l + 3 sin 2 9)i] as + 2 cos 3 9n ss ] sign(6^)sign (^^j 



«c = 



tine (in^ 



+ tanA 



4?? CTS (?? CTCT -Vss) ~ ^las + ^ sin 2 sin 20 cos 26>e fe e - ^ sin 3 6* cos 9^ h e 



dink dink 

+ | sin 9(3 cos 2 9 — l)r] ss — 5 sin 2 9 cos 9r] as ^ sign(6 v ,)sign 



V 



'e b e 



+ tan 2 A 



4r/ 2 s - Arj as y/e^e sin 6 9 - sin 9 ebe \ , 



6 , 



(3.2) 



where 

dn^ 



dink 



dns 
din i 



6=0 



=2(-7 + 4 cos 2A - cos 4A)e 2 + 2(4-3 cos 2A + cos 4A)er/ CT(T + 4(2 sin 2A - sin 4A)er] as 

+ 2(2 + cos 2A - cos 4A)er/ ss + sin 2 2A(r/ 2 (T + rf ss ) + 4 cos 2 2A?? 2 S - (1 - cos 4A)r/ (7(T r ?s 
+ 2 sin 4A (r^r^ - VcTS Vss) - (1 - cos2A)e CT(JCT - 2sin2A£ 2 CTS - (1 + cos2A)£ 2 ss , 



= - 8e 2 + 4e(r/ CT(T + rj ss ) + 4r/ 2 s - 2& ss . 



6=0 



(3.3) 



All slow-roll parameters are evaluated at horizon crossing. We note that the running 
spectral index of isocurvature perturbation, as, is independent of A, which means that the 
as is determined at horizon crossing and does not change thereafter. 
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a,s 



H T ss A 



s x %>x ^xx 

@ Vera Vers Vss 



e 2 



i 



XXX 
2 



e 2 e 2 

t2 p2 



Table 1: Full set of model parameters to second order, listed by physical fields (y,x) or component 
fields (<j,s). H, Tss and CLs are n °t required for the spectral running and can be disregarded. 



4. Non-Gaussianity During Inflation 

Another potential discriminator between inflationary models is the level of non-Gaussianity 
produced during inflation. Current observations limit the nonlinearity parameter, |/nl| < 
100 [30]. A perfect CMB experiment cannot hope to detect |/nl| < 3 [31]. Using the SN 
formalism [32, 33, 34, 35, 36], the non-linearity parameter is given by [37] 

_6 A IJK N tI N t jN, K g IM G KN N I N K N M N , n 

5 /NL (N I N J g^)^ i kf + (NjNjguy ' { - } 

fi (3) 

which is valid in the slow-roll limit. The first term, which we call — §/nl, can be written 
as [38], 

-f/S = ^(l + /). ( 4 - 2 ) 

where r = Vt/Vq is tensor-to-scalar ratio and / is a function of the momentum triangle 
with the range of values < / < | [16]. In Appendix C, we have shown Eqn.(4.2) is valid 
for non-canonical kinetic terms. From [28], the tensor-to-scalar ratio is calculated to be 
r < 16e and therefore we can approximate 



6 f (3) 

" g J NL 



< 



(4.3) 



which is certainly too small to be observable in CMB experiments. Also see [39]. Due to 
this, we will concentrate on the second term, only in the next sections and consider 
two cases: separable potentials of product and sum. 

4.1 Product Potential W(<p,x) = U{ip)V(x) 

We first consider the case of a separable potential by product, for which we derive the 
analytic formula for the non-linear parameter /nl- To use Eq. (4.1), we need to know the 
dependence of the number of e-foldings on the fields <p and %. Following the calculations 
in [6], we can obtain the derivatives of N by and %*. 
First we find the number of e-foldings, 

'-wj,h v = 'TPS (4 ' 4) 

where superscript (or subscript) e and * denotes the values evaluated at the end of inflation 
and horizon crossing respectively. From the equations of motion, in the slow-roll regime, 
it is possible to find a constant of motion along the trajectory, as in [6]: 

Cl = -— /V 2 ^)^ + — (4 5) 

1 _ Ml J Ml J V x - l4 -" j 
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Using this constant of motion, in the slow-roll regime we find the first derivatives of N(t e , t* 
with respect to the fields, 



dN sign(U/U^ 



dN _ sign{V/V x ) 



1 _ _X g 2b e -26» 



X c 2 fe e — b, 



(4.6) 



From the first derivatives we can find the second derivatives, 
Mp- 



r2 d 2 N 



Mp 



Mi 



where 



d 2 N 
dx 2 * 
d 2 N 
d(p*dx* 

Vss 

Ap 



1 _Hk + Lign^sign e l e l e ^-2K _ l e 4 be -4^ A 



,26, 



1 - 



XX 



fx e 2fc e -26» _ J_ g 4b e -46 

ze x/ e e x 



sign — f- sign — y 



U 



( e x r l<p'p + e v r ?xx) 



Vv 



1 



^ / x/44 



4b e- 3b *Ap, 



(4.7) 



e e 
V X 



1 



(ei) 2 



(4.8) 
(4.9) 



We have followed the notation of [38] as much as possible, for reasons that become apparent 
in the next section. 

From Eqs. (4.6) and (4.7) we find the second term of Eq. (4.1), 



_^f( 4 ) 

g J NL 



2g— 2fe e +26* 



u 3 a 3 



1 



2e; 



+ 



+-sign(6^)sign 



vu 2 a 2 

IT 



(e*) 2 v 6 v 



x _!?xx 
2e * 



V e* 



2b e —2b* 



A P 



with the definitions 



u = 



v = 



a = e 



-26 e +2£>* 



1 + JL fl _ e 2fc e -26«\ 



,(4.10) 



(4.11) 



When &(<£>) = 0, then a = 1 and ej, = 0, which means that the symmetry between 92 and x 
in Eq.(4.10) is restored. 

4.2 Sum Potential W(p,x) = U(<p) + V(x) 

The second case we consider consists of separable potential models by sum. Similar cases 
with canonical kinetic terms were previously studied in [38], where it was shown that they 
do not generate significant non-Gaussianity. We will now derive the general formula for 
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/nl i n the presence of non-canonical kinetic terms, closely following the derivation found 
in the above-mentioned paper. 

The total number of e-foldings along a trajectory in slow-roll regime is given by 

We note that in principle <p in the integration of x can be re- written in terms of ip*, \* 
and x along the trajectory using equations of motion. For sum potentials, ip and x have 
the relation along the trajectory through 

"X i 



rip e -2%>) rx 

/ Ij — ** = / 

Jip* U <P Jx* 



x^x 



dx- 



(4.13) 



As before, the integral of motion along the trajectory leads to a constant of motion 



= -M 2 P j , 



This enables to derive the first derivatives of N 



M P 
M P 



ON 

dip* 

dN 



sign 



sign 



Up 

w 

Yx 
w 



1 



(u* + z e )-g, 



(v* - z e ) - n, 



using the definitions, 



0((p*,x* 



M P 



Vy Qlf* 



In the same way the second derivatives can be derived: 



Mi 



Mi 



Ml 



d 2 N 
dip 2 * 
d 2 N 
dx* 
d 2 N 
dx*dip* 

d 2 N 
dip*dx* 



1 - 



2e* 



W* 



+ sign -f- 



U,, 



Mr 



dZ e 



W J W*J2e%dip" 



-M P 



do_ 

dip* 



r? xx V* - Z e 



2e* 



u <f 

sign I — - 



W* 
M P dZ e 



sign — - 



V y \ Mpe~ h * dZ e 



W J W\/2Kdx* 



M P 



W I W\/2e*dx* 

' V A. 

dG 



M P 



&H 

dx* 



dx*' 



— J>* 



-sign (&„,) sign 



w 



'el V* 



e~ W* 



— sign 



Mr 



dZ £ 



W J W*J2lldip* 



(4.14) 



(4.15) 

(4.16) 
(4.17) 

(4.18) 



- M P 



(4.19) 

on 



dip* 



The order of the second derivatives, 



d 2 N 



d 2 N 



Q — Q dip , and the last two equations emphasise 
that Q and TL are not independent quantities. In fact, they are related by 



M P 



do_ 

dx* 



M P 



&H 

dip* 



sign(6 (/P )sign -f 



J>* y* 



W 2 W*' 



(4.20) 
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A similar expression relating J^- to J^- and can be obtained, but it is irrelevant for 
our purposes and is therefore not given. From the definition of Z e in Eq. (4.16), we can 
calculate 



^^w)^wrw P w ^ As+Bs+c ^ (421) 



sign [w ) r-w, = ~i^ w * [As + Bs] ' (422) 

where we define 

- -wrf i 1 - f - (^) v^i) e " > ^ 4 - 23 ) 

6, ee ^sign(6,)sign J ^^e 2 *- 26 - , (4.24) 

^ - -^ign(6^)sign (4.25) 
From the results of the first and second derivatives, it is straightforward to calculate 



6 f (4) = 2 

C-'NL / 2 2 2 2 \ 2 

4 



/ M 2 a^ , t>V : 

\ e* ' e* 

\ V X 



4 V 2e £ d f* e <p 



where we have defined 

U* + Z e V* - z e ti . 

sign (^) (|) ^^W. (4.28) 

When 6(99) = 0, this calculation of simplifies, since 

#s = C s = Q = H = 0, a u = a v = 1, (4.29) 
and the result is identical to that found in [38]. 

5. Examples of Scalar- Tensor Theories with Product Potentials 

At first, we consider a separable potential by product with exponents in the ip field, as in 
the case of a massless dilaton, ip, 

W(<p, X ) = U(<p)V( X ) = e 4c ^V( X ). (5.1) 
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Model 


la lb 


2a 2b 


v{x) 

Vxx 
^xxx 

6 


-3-2- -3-2- 

^ M P > J M P 

|x 4 k m xX 2 

2 fc X e X 

14 o 




fx 4 i«4* 2 

2 e X e X 

14 

-16/3 -16/3 



Table 2: Four Jordan-Brans-Dicke models considered in the text and their slow-roll parameters for 
a product potential, W = e 4b ^V(x)- In these specific models, 5 parameters are reduced further to 
only 2 and 3 independent parameters for model 1 and 2 respectively. e v and e x are directly related 
to e through 6. The three remaining parameters {f] xx , £, xxx j £&) are shown to be also functions of 
e through e v and e x . 



When considering the physical fields (<p,x)> only H and T$s can be discarded. Therefore, 
the set of parameters consists of A and 11 slow- roll parameters as shown in Table 1. 
However, with the potential chosen above, the set reduces to only 5 independent slow-roll 
parameters, which are useful for the running spectral indices: 

T 2 2 _ _ Mp fV x \ nb 



= e x = ^f j e-» (5.2) 

Vxx = Ml i^f) e~- H xx = M P (^) e"-, (5.3) 

£ c = 8M|c w , (5.4) 
since the remaining 8 are not independent: 

€c = tip i (5-<5) 

V w = + 2e^ , 7]^ = 2sign(c^)sign V e <p e x = e sin26>, (5.6) 

Ctptptp = ^ e <p + 3£, c £tp + £,cc\[tip i ippxx = ^ e <p e x e x^ c ' (5-7) 

£pxxx = 2sign(c^)sign (^f) ^^i] xx . (5.8) 
5.1 Jordan-Brans-Dicke theory 

First we apply our results to the Jordan-Brans-Dicke (JBD) theory with quadratic and 
quartic potentials. In this theory, 

£%)=e 4 ^), bfr) = -l3JZ-, (5.9) 



where we assume 3 is a positive constant. For this potential choice, = e& = rj tp(f /2 = 83 2 
and the slow-roll ends when e e = e e x + = 1. 
Two potentials are chosen for %: 
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. Model la V( X ) = %X* 



• Model lb V( X ) = \m\x 2 



where A is a dimensionless parameter and m x is the mass of the x held. The slow-roll 
parameters for both quartic and quadratic potentials for the x~ ne ld are shown in the first 
two columns of Table 2. As seen in the table, the choice of potentials reduces the 5 
independent slow-roll parameters to only two: e v and e x . Equivalently, since e v = ecos 2 # 
and e x = esin 2 6, we can use the parameters e and 9. Hence, all the first and second order 
slow-roll parameters are proportional to e and e 2 respectively which make the primordial 
spectral indices proportional to e 2 in the lowest order. 

One final parameter is required, tan A, to describe the evolution after inflation. Then 
e 1 "d'lnfc S> are J US ^ a mnc ti° n °f $ an d A. We note that 9 is calculated at horizon crossing 
and depends on (3 and the initial values of <p and whereas A represents the evolution 
after horizon crossing and depends on the late time evolution of the Universe. 

In Figure 1 we show the running spectral indices and «c for Model la. The result 
for Model lb is almost identical and is not shown. As said before, as is independent of 
A, thus we plot the as dependence on 9 separately in Figure 3. We plot in the range of 
< 9 < it since the potential has the reflection symmetry x ~~ * ~~ X which corresponds to 
the change of 9 to —9. 

For the correlated spectral running, ac, there is an divergence at A ^ |. This 
divergence is not observable. From Eq. (2.23), it can be seen that, for A = | (T^g = 0), 
there is no evolution in and the amplitude of the correlation spectrum is zero. It is 
therefore impossible to define a spectral index or running at this point. As A — ► ^, T^s is 
non-zero, but very small. In this limit, the correlation amplitude would be also small and 
therefore unobservable. 

As an aside, it can be noted that, for the specific case of Model la, B (defined in 
Eq. (2.19)) is identically zero, due to a cancelling of the slow-roll parameters. This means 
that from Eq. (2.18) the curvature perturbation remains constant after horizon crossing 
(C = 0) for this JBD model with quartic potential. 

It is possible to consider the case V{x) = ^X n - in this JBD theory, and using the fact 
that at the end of inflation e e = e^, + t e x = 1, we can analytically solve iV in the slow-roll 
limit, 



1 + 4 ff -2/WAfp 2 

L+ nM 2 p 6 X * 



1 



1 + (2n - 8)/3 2 



(5.10) 



We find that (3 < 0.05 is required to obtain enough e-foldings for tp*,x* < 50Mp. From 
this fact, 8(3 2 N < 1 and the limit of Eq. (5.10) gives e* ~ < 0.004n, where we have 
taken iV k, 60. In this limit of (3, we can approximate Eq. (5.10) further: 

N m — — — — (5.11) 

2n 4 ' n y 1 

For Models la and lb, we see from Figure 1 that 

-10(6*) 2 < a c < 0. (5.12) 
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(n-1) / 6 



< TT 




< 7T/2 



7T/2 7T 



<a tt/2 



a r / e 




-3.0 -2.5 -2.0 



-1.0 



Figure 1: The spectral index ({n^ — l)/e) and running of ri( and nc in terms of e for model 1-a. 



Hence \a^\ < 1.6n 2 x 10 4 which is quite compatible with observations [14]. 

From the slow-roll equations of motion, Eq. (2.16), we find the solution of ip 



r ~ v" + 3 //-si oil ( -jj- 



'2e v N, 



(5.13) 



which gives b e 
(we take e e 



e , e 



-4(3 2 N. Using this we can express Eq. (4.10) at the end of inflation 

1) in terms of slow-roll parameters at horizon crossing. In the limit 

8P 2 N <C 1, we find a as 1 + N, u e e and v « 1 — e e . This leads to the simplified form of 
f (4) 

JNLi 



--f (4) ~2e* 



"XV 



(5.14) 



It can be noted that in this small (3 limit, e* < e* and the last two terms in this equation 
are negligible compared to the first two terms. 

Using Eq. (5.11), can be written in terms of the required e-folding number in the 
small (3 limit (required for enough inflation) 

P P* 



6 ,(4) _ _•_ 
5 Tnl ~2N NM P ' 



(5.15) 



It is clear that, in order to achieve enough inflation, the parameters are required to be 

2]y, which is unobservable. 



small and hence — §/jyx 



5.2 Brans-Dicke Type Models with Quadratic Exponent 

In order to observe the effect of higher powers in b((p), we consider a quadratic function, 
such that ^ 0. To this end, we consider a Brans-Dicke type model with 

l%) = e 4fe ^, %) = -/V7Mj>. 

With this choice, e v = 32 / #V/M| and r/ w = -8/3 + 64/?V/^J>- 
Once again, two example potentials are chosen for x- 



(5.16) 
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-3.8 -3.1 -2.4 -1.7 -1.0 -30 -25 -21 -16 -12 -8 -3 5 

Figure 2: Spectral index (n^ — 1) and running of tiq and nc for model 2-a. We used (3 = O.le here. 




Figure 3: Running of ns for Models 1 and 2. as is independent of A. 

. Model 2a V( X ) = \x 4 

• Model 2b V(x) = \rn 2 x \ 2 

where A and m x have the same definitions as before. 

The slow-roll parameters for these potentials are shown in the last two columns of 
Table 2. Once again, most of the slow-roll parameters can be written in terms of the single 
parameter, e 2 . The exception is £&, which is related to the coefficient of the coupling term, 
(5. In order to analyse the expected levels of running from these cases, it is necessary 
to relate j3 to e and we take the assumption that (5 = O.le. For Model 2a, the numeric 
calculations of the spectra for this special case (in relation to e, 9 and A) are shown in 
Figures 2 and 3. Model 2b is almost identical to Model 2a and the numerical results are 
not shown. 

Models 2a and 2b have a reflection symmetry about <p (0 — > —9) and/or x (9 ^ tt — 9). 
However since the sign of tan A is determined by the sign of a in Eq. (2.20) and a is 
proportional to sin 29 for the potential in Eq. (5.1), The reflection of ip or x changes also 
the sign of tan A before the end of inflation. For this reason, we see the symmetry in 
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Figure 2 around the point {6 = 7r/2,A = tt/2). However after inflation, A is no longer 
dependent on the fields cp and x, thus this symmetry is not necessarily valid any more. 
In this model, the solutions for <p and e^, are given by 

= 4 = W\l/Ml = e;e l ^ N . (5.17) 

Due to the quadratic coupling, the potential is steeper than that of traditional Brans-Dicke 
models, so that a much smaller 3 is required to give enough e-folding number. In this case, 
if we assume 16/3 N <C 1, then we obtain 



1 + 8 le-W- x l 
n 



1 In 



16/3 



ABn 

1 + 



1 - 32/3 V 2 J ' 



(5.18) 



In order to obtain enough inflation, for ip*,x* < 50Mp, then we require (5 < 0.0005. In 
this very small limit, we can approximate Eq. (5.18): 



2n 4 +P \ 2 n n 2 1 " 1 j 



As in Model la, it is possible to estimate the non-Gaussianity for small 3. Again, we 
take e e = 1, so that a, u and v are given as before. In the limit 16/3/V <C 1, we find 

-\ fSl - K - <x + K^ N + 1) - 4 - \ (4) 2 W + x ) 3 - < 5 - 2 °) 

The final term is directly due to b w / 0. Using the approximation in Eq. (5.19), we find 

_V 4) ~ — + (5 21) 

5 1nl- 2N + NM 2 p 

As in Model 1, in order to achieve enough inflation, the parameters are required to be small 
and hence -f/^ ~ j N - 

6. Examples of Scalar- Tensor Theories with Sum Potentials 

The second models we consider are those with sum potentials. If we define the potential 



W(<p, X ) = U(ip) + V{x) = \m 2 y + l -m 2 xX \ (6-1) 



then the slow-roll parameters simplify greatly and are given by 

2 2 

r? w = M 2 ^, Vxx = M^e- 2b , r, m = 0, (6.3) 
£vw ~~ ^ffXX ~ £<pxxx ~ ^xxx ~ ^' 
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Figure 4: Spectral index (ri( — 1) and running of and «c in terms of e for model 3. Values of 
(3 = 0.1 y/e and re~ 2b = 2 have been assumed. 




Figure 5: Running of ns for Models 3 and 4. as is independent of A. 

It is no longer possible to directly relate the parameters to e and 9. However, if we 
represent the ratio of masses by 



(6.5) 



then the slow-roll parameters can be re-written in terms of e, (3 and re only: 
2m, , sin2 W 



f? w = e I cos (9) + 



re 



-26 



£6 = 64 



r/ xx = e (W(#) + cos 2 (#)r e - 2fc 
/3 2 cos 2 (fl) 



(6.6) 



(6.7) 



For this sum potential, the remaining choice is in the function of b(ip) and we consider 
the two obvious cases: 



Model 3 %) = 
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• Model 4 b(<p) = -ff£ 

We can fix (3 and the mass ratio, in order to estimate the spectral index and running. 
We take re~ 2b = 2 and, for Models 3 and 4 respectively, we assume [5 = O.ly^e and f3 = O.le. 
The results of the spectrum for Model 3 are shown in Figures 4 and 5. The results for 
Model 4 are very similar to those of Model 3 and are therefore not shown. 

7. Numerical Analysis of Non-Gaussianity 

In order to verify our analytical results of the previous section, we numerically solved the 
background equations of motion, Eq. (2.1) and the Friedmann equation, Eq. (2.2), not 
assuming slow-roll. Multiple trajectories were considered, with initial values, ipi and Xii 
given by a (p-\ grid and <fi = Xi = 0. The end of inflation is taken to be the time at which 
the slow-roll parameter e = + e x becomes unity and the total number of e-foldings until 
this point is given by Nf. 

At horizon crossing, the fields have reached values of 99* and x* an d the remaining 
number of e-folds is denoted by N. The calculation of non-Gaussianity requires the gradient 
of the number of e-foldings between horizon crossing and the end of inflation, N with respect 
to each of the fields at crossing, 93* and %*. 

It is useful to notice, quite trivially, that if for a given trajectory, </?j = ip* and Xi = X*> 
then Nf = N. It is therefore possible to calculate N,i, N Y ij etc by differentiating the grid 
of Nf with respect to ipi and Xi- 

Each of the cases considered above (Models 1-4) were modelled numerically. The 
analytical slow-roll approximations given in Eq. (5.14) and (5.20) were seen to agree within 
5%. For Models 3 and 4, no analytical approximation was found. For these sum potentials, 
the numerical level of non-Gaussianity was calculated from the second term in Eq. (4.1) 
and the results are shown in Figure 6. 

In previous work, [38], for theories with canonical kinetic terms and a sum potential, 
/nl ~ l/N. We find numerically that the /nl is of that order of magnitude, even with 
non-canonical couplings. Therefore we conclude that (3 has little impact on the level of 
nonlinearity. 

8. Conclusion 

The focus of this paper was slow-roll inflation in theories with two scalar fields, coupled 
through a potential as well as their kinetic terms. We have derived the general formulae for 
the running of the spectral indices for adiabatic and entropic perturbations as well as for 
the cross-correlation spectrum. In addition, using the (5iV-formalism and specialising to the 
product potential {W(ip, x) = U ( l f)V(x)) and the sum potential {W(ip, x) = U(ip) + V(x)) 
we have derived the general expressions for the nonlinearity parameter /nl during inflation. 

We have calculated the spectral indices and runnings for several example models. The 
results are within present observational ranges. 

One of the specific examples we have considered is the Jordan-Brans-Dicke theory with 
quadratic and quartic potentials. In this case, the slow-roll parameter for the field 99 is 
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Figure 6: Level of non-Gaussianity for Models 3 (left) and 4 (right) calculated from the second 
term of Eq. (4.1). We have assumed m v = m x = 5x lO~ 5 Mp. (3 = 0.1. 

given by the (constant) coupling parameter j3, denned in Eq. (5.9). The slow-roll condition 
for the field ip implies that (3 itself has to be small. In this case, an analytical formula for 
/nl during inflation is given by Eq. (5.14) and is of order 1/N, where N is the total number 
of e-foldings. Assuming to be positive, the effect of the coupling is to enhance /nl by 
an additional term jffif- ■ When the coupling is quadratic in the yj-field (as in Eq. (5.16)), 

the nonlinearity is enhanced by a term j^jjt ■ Neither of these terms are significant and in 
the slow-roll limit, /nl is not much enhanced by the coupling. 

We also considered models with a sum of two quadratic potentials. Due to the com- 
plexity of the analytic formula for /nl> it is difficult to find an approximate equation. We 
therefore provide numerical results for simple models and showed that /nl "C 1 and is 
unobservable for the slow-roll limit. 

To conclude, in order to distinguish between inflationary models, future observations 
will search for deviations from standard one-field inflation. In order to do so one has to go 
beyond the spectral index as an observable. Future experiments, such as PLANCK, will 
tighten the constraints on the spectral running and nonlinearity. The results of this paper 
can be used to calculate these parameters in other multi-field models with non-canonical 
kinetic terms, in order to constrain the parameter ranges. 
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A. Slow-Roll Parameters 

The slow-roll parameters we shall use are first-order 
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ln(l + T£ s ) = 2 sin A cos A(-a, - 5*T (S ) 



^T SS = -K (B.5) 

d . _ -a* - 5*T CS . 
■lni C 5 = 



C. Calculation of — for Non-Canonical Terms 

The first term of Eq. (4.1), /®, can be written by 



6 (3) _ A IJK N T NjN K 
V NL (NjNjg")^^ 



— f£L= ,Cr *r /^V,^ (CI) 



where summations over /, J, K are implied and .A 7 is given by [37] 



P perms 



(C.2) 



with momentum dependent M123 = M(ki, k 2 , k 3 ). We can perform the sum and permu- 
tations in the numerator of Eq. (C.l): 

E A IJK NiNjN K = E NjNjNkjL E ^^^^123- (C.3) 

UK UK p perms 

Following from this, the numerator in Eq. (C.l) can be reduced to 

-4^2 (E^^^J E (C4) 
P \ U / perms 



where we have used Yli Njfii = —H. We define 



y 



f = f( kl ,k 2 , k 3 ) = -1 - ^ (C.5) 
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Finally, if we consider the power spectra of curvature perturbations and gravitational waves, 
we obtain 

-f/S = ^(i + /). (c- 7 ) 

where r is scalar-to-tensor ratio, r = V g /V(. 



-22- 



References 



E. J. Copeland, J. Gray and A. Lukas, Moving five-branes in low-energy heterotic M-theory, 
Phys. Rev. D 64 (2001) 126003 [arXiv:hep-th/0106285]. 

P. Brax, C. van de Bruck, A. C. Davis and C. S. Rhodes, Cosmological evolution of brane 
world moduli, Phys. Rev. D 67 (2003) 023512 [arXiv:hep-th/0209158]. 

L. Cotta-Ramusino and D. Wands, Low-energy effective theory for a Randall- Sundrum 
scenario with a moving bulk brane, arXiv:hcp-th/0609092. 

D. Wands, N. Bartolo, S. Matarrese and A. Riotto, An observational test of two-field 
inflation, Phys. Rev. D 66 (2002) 043520 [arXiv:astro-ph/0205253]. 

J. Garcia-Bellido and D. Wands, Constraints from inflation on scalar - tensor gravity 
theories, Phys. Rev. D 52 (1995) 6739 [arXiv:gr-qc/9506050]. 

J. Garcia-Bellido and D. Wands, Metric perturbations in two-field inflation, Phys. Rev. D 53 
(1996) 5437 [arXiv:astro-ph/9511029]. 

V. F. Mukhanov and P. J. Steinhardt, Density perturbations in multifield inflationary models, 
Phys. Lett. B 422 (1998) 52 [arXiv:astro-ph/9710038]. 

S. Tsujikawa, D. Parkinson and B. A. Bassett, Correlation- consistency cartography of the 
double inflation landscape, Phys. Rev. D 67 (2003) 083516 [arXiv:astro-ph/0210322]. 

A. A. Starobinsky, S. Tsujikawa and J. Yokoyama, Cosmological perturbations from 
multi-field inflation in generalized Einstein theories, Nucl. Phys. B 610 (2001) 383 
[arXiv:astro-ph/0107555] . 

P. R. Ashcroft, C. van de Bruck and A. C. Davis, Boundary inflation in the moduli space 
approximation, Phys. Rev. D 69 (2004) 063519 [arXiv:astro-ph/0310643]. 

C. Ringeval, P. Brax, C. van de Bruck and A. C. Davis, Boundary inflation and the WMAP 
data, Phys. Rev. D 73 (2006) 064035 [arXiv:astro-ph/0509727]. 

S. Tsujikawa and B. Gumjudpai, Density perturbations in generalized Einstein scenarios and 
constraints on nonminimal couplings from the Cosmic Microwave Background, Phys. Rev. D 
69 (2004) 123523 [arXiv:astro-ph/0402185]. 

C. T. Byrnes and D. Wands, Curvature and isocurvature perturbations from two-field 
inflation in a slow-roll expansion, Phys. Rev. D 74 (2006) 043529 [arXiv:astro-ph/0605679]. 

D. N. Spergel et al, Wilkinson Microwave Anisotropy Probe (WMAP) three year results: 
Implications for cosmology, arXiv:astro-ph/0603449. 

N. Bartolo, E. Komatsu, S. Matarrese and A. Riotto, Non-Gaussianity from inflation: 
Theory and observations, Phys. Rept. 402 (2004) 103 [arXiv:astro-ph/0406398]. 

J. M. Maldacena, Non-Gaussian features of primordial fluctuations in single field inflationary 
models, JHEP 0305 (2003) 013 [arXiv:astro-ph/0210603]. 

G.I. Rigopoulos, E.P.S. Shellard and B.J.W. van Tent, Quantitative bispectra from multifield 
inflation, arXiv: astro-ph/0511041. 

T. Battcfcld and R. Easther, Non-gaussianities in multi-field inflation, 
arXiv:astro-ph/0610296. 



- 23 - 



[19] X. Chen, R. Easther and E. A. Lim, Large non-Gaussianities in single field inflation, 
arXiv:astro-ph/0611645. 

[20] L. Alabidi, Non-gaussianity for a two component hybrid model of inflation, JCAP 0610 
(2006) 015 [arXiv:astro-ph/0604611]. 

[21] N. Bartolo, S. Matarrese and A. Riotto, Non-Gaussianity from inflation, Phys. Rev. D 65 
(2002) 103505 [arXiv:hep-ph/01 12261]. 

[22] F. Bernardeau and J. P. Uzan, Non-Gaussianity in multi-field inflation, Phys. Rev. D 66, 
103506 (2002) [arXiv:hcp-ph/0207295]. 

[23] K. Enqvist, A. Jokincn, A. Mazumdar, T. Multamaki and A. Vaihkoncn, "Non-Gaussianity 
from Preheating," Phys. Rev. Lett. 94 (2005) 161301 [arXiv:astro-ph/0411394]. 

[24] A. Jokinen and A. Mazumdar, "Very Large Primordial Non-Gaussianity from multi-field: 
Application to Massless Preheating," JCAP 0604 (2006) 003 [arXiv:astro-ph/0512368]. 

[25] X. Chen, M. X. Huang, S. Kachru and G. Shiu, Observational signatures and 
non-Gaussianities of general single field inflation, arXiv:hep-th/0605045. 

[26] N. Barnaby and J. M. Cline, Nongaussian and nonscale-invariant perturbations from 
tachyonic preheating in hybrid inflation, Phys. Rev. D 73, 106012 (2006) 
[arXiv:astro-ph/0601481]; N. Barnaby and J. M. Cline, Nongaussianity from tachyonic 
preheating in hybrid inflation, arXiv:astro-ph/0611750. 

[27] F. Di Marco, F. Finelli and R. Brandenberger, Adiabatic and Isocurvature Perturbations for 
Multifield Generalized Einstein Models, Phys. Rev. D 67 (2003) 063512 
[arXiv:astro-ph/0211276]. 

[28] F. Di Marco and F. Finelli, Slow-roll inflation for generalized two-field Lagrangians, Phys. 
Rev. D 71 (2005) 123502 [arXiv:astro-ph/0505198]. 

[29] C. Gordon, D. Wands, B. A. Bassett and R. Maartens, Adiabatic and entropy perturbations 
from inflation, Phys. Rev. D 63 (2001) 023506 [arXiv:astro-ph/0009131]. 

[30] P. Creminelli, L. Senatore, M. Zaldarriaga and M. Tegmark, Limits on /nl parameters from 
WMAP Syr data, arXiv:astro-ph/0610600. 

[31] E. Komatsu and D. N. Spcrgel, Acoustic signatures in the primary microwave background 
bispectrum, Phys. Rev. D 63, 063002 (2001) [arXiv:astro-ph/0005036]. 

[32] A. A. Starobinsky, Multicomponent De Sitter (Inflationary) Stages And The Generation Of 
Perturbations, JETP Lett. 42, 152 (1985) [Pisma Zh. Eksp. Teor. Fiz. 42, 124 (1985)]. 

[33] M. Sasaki and E. D. Stewart, A General analytic formula for the spectral index of the density 
perturbations produced during inflation, Prog. Theor. Phys. 95 (1996) 71 
[arXiv:astro-ph/9507001] . 

[34] M. Sasaki and T. Tanaka, Super-horizon scale dynamics of multi-scalar inflation, Prog. 
Theor. Phys. 99, 763 (1998) [arXiv:gr-qc/9801017]. 

[35] D. H. Lyth, K. A. Malik and M. Sasaki, A general proof of the conservation of the curvature 
perturbation, JCAP 0505, 004 (2005) [arXiv:astro-ph/0411220]. 

[36] D. H. Lyth and Y. Rodriguez, The inflationary prediction for primordial non-gaussianity, 
Phys. Rev. Lett. 95 (2005) 121302 [arXiv:astro-ph/0504045]. 



- 24 - 



[37] D. Seery and J. E. Lidsey, Primordial non-gaussianities from multiple-field inflation, JCAP 
0509 (2005) Oil [arXiv:astro-ph/0506056]. 

[38] F. Vernizzi and D. Wands, Non-Gaussianities in two-field inflation, JCAP 0605 (2006) 019 
[arXiv:astro-ph/0603799] . 

[39] I. Zaballa, Y. Rodriguez and D. H. Lyth, "Higher order contributions to the primordial 
non-gaussianity," JCAP 0606 (2006) 013 [arXiv:astro-ph/0603534]. 



- 25 - 



